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Abstract: The large N limit of the anomalous dimensions of operators in A/" = 4 super 
Yang-Mills theory described by restricted Schur polynomials, are studied. We focus on 
operators labeled by Young diagrams that have two columns (both long) so that the 
classical dimension of these operators is 0{N). At large these two column operators 
mix with each other but are decoupled from operators with n ^ 2 columns. The planar 
approximation does not capture the large N dynamics. For operators built with 2, 3 
or 4 impurities the dilatation operator is explicitly evaluated. In all three cases, in a 
certain limit, the dilatation operator is a lattice version of a second derivative, with 
the lattice emerging from the Young diagram itself. The one loop dilatation operator 
is diagonalized numerically. All eigenvalues are an integer multiple of Sqym there 
are interesting degeneracies in the spectrum. The spectrum we obtain for the one loop 
anomalous dimension operator is reproduced by a collection of harmonic oscillators. 
This equivalence to harmonic oscillators generalizes giant graviton results known for 
the BPS sector and further implies that the Hamiltonian defined by the one loop large 
dilatation operator is integrable. This is an example of an integrable dilatation 
operator, obtained by summing both planar and non-planar diagrams. 
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1. Introduction 

In the last few years interesting new progress has been made in the study of the dy- 
namics of multimatrix models. Starting from the remarkable observation[l] that the 
Schur polynomials are a complete basis of gauge invariant operators which diagonalize 
the two point function of the free {gy^ = 0) super Yang-Mills theory, similar bases 
have been found for multimatrix models[2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. For these 
bases, the two point function is diagonal and known exactly as a function of N (but of 



- 1 - 



course, at Qym = 0)- The fact that the dependence is known exactly suggests that 
these results will be useful for going beyond the planar approximation. 

When would this be needed? Often "the planar hmit" and "the large N hmit" 
are taken as synonyms. This is not, in general, accurate. For example, if we imagine 
computing the two point correlator of an operator with a bare dimension A of most 
A ~ J with ^ <^ 1, then summing the planar diagrams will capture the large N 
limit. For operators with a dimension larger than this, exploding combinatoric factors 
overpower the nonplanar (-^) suppression and the planar approximation is completely 
ineffective [13]. In this scenario, to get the correct large N limit, it is necessary to sum 
a lot more than just the planar diagrams. On general grounds we expect the large N 
limit to be simpler than the full theory[14]. The planar diagrams are a small subset 
of all possible diagrams, so that it is quite natural to expect that summing only the 
planar diagrams will give a much simpler problem. Why should the large N limit be 
simple when one needs to sum much more than just the planar diagrams? The answer 
to this question will probably not be general, but rather will depend on the specific 
dynamical problem considered and must be answered case by case. A very pedestrian 
approach is simply to compute the large limit and then to look for simplifications. 
This has been accomplished[15, 16, 17] in a number of interesting examples including 
LLM geometries [18, 19] and the near horizon geometry of a bound state of giant 
gravitons[20]. The results are remarkably simple. Indeed, as an example, for ^ BPS- 
correlators in the presence of M giant gravitons with M of order N, [15, 16, 17] showed 
that the usual expansion is replaced by a expansion. Further, if one expands 
the exact correlators (which because they are | BPS do not depend on but only 
on or + M) the expansion coefficients for correlators in the background of M 
giants are exactly the same as the expansion coefficients for correlators with no giants 
present! This remarkably simple result was confirmed holographically[21] by matching 
to graviton dynamics in the LLM geometries using the formaUsm of [22]. For near-EPS 
operators corresponding to BMN loops[23] it was argued in [24, 25, 26] that the usual 
't Hooft couphng Qym^ replaced by the effective 't Hooft couphng Qyui^ + 
For additional interesting related studies see[20, 27]. 

In this article we will consider the problem of computing the anomalous dimension 
of an operator with a bare dimension of order A^. To answer this question, we need to 
go well beyond the planar limit; we find the methods and approach of[2, 3, 4, 5, 6, 7, 
8, 9, 10, 11] once again surprisingly effective. The results again exhibit a remarkable 
simplicity - the spectrum of anomalous dimension can be matched to the spectrum of 
a set of oscillators! We see once again that the large A'" hmit is indeed a simple hmit. 

The operators we consider, restricted Schur polynomials, will be built using O(A^) 
Zs and 3 or 4 "impurities" {Ys), where Z and Y are complex adjoint scalars ol N — ^ 
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super Yang Mills theory. The case of operators with two impurities was studied in[30]. 
The dilatation operator when acting on a restricted Schur polynomial, produces terms 
that have a combination ZY — YZ appearing. In [30], the techniques of [4, 6] were 
used to separate the Z and the Y and then write the results as a linear combination 
of restricted Schur polynomials. This method is very cumbersome as it involves the 
inversion of a matrix. This must be done analytically so its tedious (for the case of two 
impurities one must invert a 6 x 6 matrix). More than two impurities was effectively 
out of reach. In this article we develop a new formula (in section 2) which avoids 
this matrix inversion. This allows us to handle the cases of three and four impurities 
without much trouble. The resulting formulas for the dilatation operator are quite 
lengthy (see Appendix A), but their spectrum is surprisingly simple. 

Our results suggest that for the class of operators considered, the Hamiltonian 
defined by the dilatation operator is integrable - it is just a set of oscillators. This is 
an example of an integrable dilatation operator, obtained by summing both planar and 
non-planar diagrams. 

The operators we consider can be mapped to giant gravitons[28] in spacetime[13, 
1, 29]. There is in fact already a known connection between the geometry of giant 
gravitons and harmonic oscillators [31, 32, 33, 34, 35, 36]. Our work differs from these 
results in at least two important ways. Firstly, we claim that the complete spectrum 
(not just the BPS spectrum!) has a connection to harmonic oscillators. Secondly, 
we have very good control over the set of operators we consider. Our operators are 
dual to a two giant system. Previous studies captured the full set of BPS states and 
consequently were not able to distinguish (for example) giant graviton plus graviton 
from excited giant graviton. Our study captures only states of the two giant system. 
Thus, it is rather natural to associate our oscillators with excitations modes of a giant 
graviton. In the same way that there are oscillators in the worldsheet theory of a string 
describing the oscillation modes of a string, we are describing the oscillators that would 
be present in a world volume description of giant gravitons describing the oscillations 
modes of a giant graviton D3 brane. 

The problem of computing anomalous dimensions for operators with a large TZ- 
charge has been considered before by a number of authors. The restricted Schur poly- 
nomials we consider in this article are built by distributing "impurities" (given by Ys), 
in an operator built mainly from Zs. If we replace the Ys by words containing 0{\/N) 
letters (which may be Z, Y or other field or derivative of a field), these words are 
naturally identified with open strings[37, 38, 39]. In this case the dilatation operator 
reproduces the dynamics of open strings ending on a giant graviton[40, 3, 4, 6]. The 
mixing of operators is highly constrained. Indeed, in [4, 6] it was shown that operators 
which mix can differ at most by moving one box around on the Young diagram labeling 
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the operator. Another interesting basis to consider is the Brauer basis[5, 11]. This 
basis is built using Brauer algebra projectors. The structure constants of the Brauer 
algebra are N dependent. There is an elegant construction of a class of BPS operators 
[41] in which the natural N dependence appearing in the definition of the operator[42] 
is reproduced by the Brauer algebra projectors [41]. Finally, another very natural ap- 
proach to the problem, is to adopt a basis that has sharp quantum numbers for the 
global symmetries of the theory [7, 9]. The action of the anomalous dimension operator 
in this sharp quantum number basis is very similar to the action in the restricted Schur 
basis: again operators which mix can differ at most by moving one box around on 
the Young diagram labeling the operator [43]. For further related interesting work see 
[44, 45]. Finally, for a rather general approach which correctly counts and constructs 
the weak coupling BPS operators see[36]. 

We now conclude this introduction with a description of what is to follow. In 
section 2 we will describe the action of the dilatation operator on restricted Schur 
polynomials. The main result of this section is formula (2.2) which gives a very explicit 
description of the action of the dilatation operator. In section 3 we will describe, in 
broad terms, the set of operators we consider and explain how the dynamics of the 
large N limit simplifies. In section 4 we explain how to construct the projectors needed 
to evaluate the action of the dilatation operator. A very detailed description of the 
specific operators we consider is given in section 5; we also describe a limit in which 
our system of two giants should be well described as a system of D3 brane giants plus 
open strings. In this limit we see that the dilatation operator reduces to a lattice 
version of the second derivative, with the Young diagram labels of the restricted Schur 
polynomials defining the lattice. In section 6 we present our numerical results and draw 
some general conclusions from them in section 7. The explicit result for the dilatation 
operator is given in Appendix A; the intcrtwiners which enter into the expression for 
the dilatation operator are described in Appendix B. 



2. Action of the Dilatation Operator 

In this section we will study the action of the one loop dilatation operator on restricted 
Schur polynomials built using two complex adjoint scalars. The main result of this 
section is the surprisingly simple result (2.2) for the action of the dilatation operator. 

We will consider the action of the one loop dilatation operator in the SU{2) 
sector [46] of A/" = 4 super Yang Mills theory 
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on the restricted Schur polynomial 

/^n,(r,S)\ -I J n\m\ ^ "^^'^ V(n+1) V(n+m) 

The labels of our restricted Schur polynomial X(R,(r,s)) are (i) i?, which is a Young 
diagram with n + m boxes or equivalently an irreducible representation of Sn+m^ (h) 
which is a Young diagram with n boxes or equivalently an irreducible representation 
of Sn and (iii) s which is a Young diagram with m boxes or equivalently an irreducible 
representation of Sm- The notation Tr(, ,,.) implies that one should only trace over the 
subspace carrying the irreducible representation^ (r, s) of Sn x inside the carrier 
space for irreducible representation R of Sn+„i- This trace is most concretely realized 
by including a projector Pfi_^(^r,s) (from the carrier space of R to the carrier space of 
(r, s)) and tracing over all of R. A simple calculation yields^ 

DxR,ir,s){Z^^Y^n= ^^f^_ J2 TV(,,,)(r^((n,n + l)^-^(n,n + l)))x 
xZ;i ...zi"-' (YZ-ZY)i" 5i"+' ...1:.^"+'" . (2.1) 

The sum over ijj runs only over permutations for which ip[n + 1) = n + 1. To perform 
this sum over ip, write the sum over Sn+m as a sum over cosets of the Sn+m-i subgroup 
obtained by keeping precisely those permutations that satisfy ■(/'(n + 1) = n + 1. The 
result follows immediately from the reduction rule for Schur polynomials (see [47] and 
appendix C of [3]) 

DXR,ir,s)= . E E ^ (^,^) (rii((^. « + l))ril' (V') 

-TR^(iP)rR((n,n + l))) Z^^ ---Zi"-^ (YZ - ZYY" . . . . 

The sum over i?' runs over all representations that can be subduced from R. Concretely, 
K runs over all Young diagrams that can be obtained from R by dropping a single box; 
crr! is the weight^ of the box that must be removed from R to obtain i?'. We will make 
use of the following notation for restricted characters 

XR,(r,s)[(y) = Tr(^,,)(^rR(a)) = Tr(^PR_^(^,,)rfl((T)) . 

Now, using the identity (bear in mind that ■0(n + 1) = n + 1) 

ZY ■■■ Z^-' {YZ-ZY)Y Y'"+^ ■ ■ ■ = Tr f f (n, n + 1) - (n, n + 1)) Z^^F®"^) 

^In general, because (r, s) can be subduced more than once, we should include a multiplicity index. 
We will not write this index explicitly in this article. 
^Our index conventions are {YZ)^. = YJZ-j,. 

^Recall that the weight of a box in row i and column j is N — i + j. 
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where 

^ ' V(n) V(n+1) V(n+m) ' 

and (this identity is proved in [48]) 

T,{t,u) 

we obtain 

^Xi?,(r,s)(^,^) = X] MR^(^r,s);T,(t,u)XT,{t,u){Z,Y) , 
T,{t,u) 

MR,(r,sy,T,(t,u) =9yM Y1 

xXT,(t,«)(V''H«>" + 1) - {n,n + l)'ip~^) . 
The sum over ip can be done by using the fundamental orthogonahty relation 

CRRfdrnm 
^ dR'dtdu{n + m 



X 



TT{{n, n + 1)), PT^(t,u) It' R' (2.2) 



The reader should consult Appendix B for a definition of the intertwiners Ir't'- This 
expression for the one loop dilatation operator is exact in N. It is one of the key results 
of this article. 

To obtain the spectrum of anomalous dimensions, we need to consider the action of 
the dilatation operator on normalized operators. The two point function for restricted 
Schur polynomials has been computed in [8] 

{XR,ir,s){Z,Y)xT,it,u){Z,Y)^) = ^R,ir,s)nt,u)fR^^^^^^^ ■ 

In this expression Jr is the product of the weights in Young diagram R and hooksR is 
the product of the hook lengths of Young diagram R. The normalized operators can 
be obtained from 



XR,ir,s,iZ,Y) = ■ 



In terms of these normalized operators 



DOR^^r,s){Z,Y) = E NR,{r,sy,T,it,u)OT,{t,u){Z,Y) 
T,{t,u) 



-6- 



R,(r,s);T,(t,u) 9ym ^ dn>dtdu{n + m)'^'fR hooksR hooks* hooks,, • 



crr> dTnm / Jt hooksy hooks,, hooks^ 



xIV ( [rR((n, n + 1)), Irt' [rT((n, n + 1)), PT-,^t,u) 



T' R' 



This last expression will be used later when we numerically study the spectrum of the 
dilatation operator. 



3. Excited Giant Graviton Bound States 

The goal of this section is to clearly define the class of operators being considered and 
to outline the approximations that can be made in the large N limit. 

In this article we will study restricted Schur polynomials labeled by Young diagrams 
with at most two columns. The number of Zs appearing is aN where 2 — a = ( <^ 1. 
The number of Ys appearing is fixed to be 0(1). These operators are dual to giant 
gravitons that wrap an in the of the AdSsxS^ background. Since the restricted 
Schur polynomials furnish a suitable basis for the two giant system, we know that 
these operators capture all excitations (BPS and non supersymmetric) of the two giant 
system. For a study of excitations of the single giant system using restricted Schur 
polynomials see [30]. For a spacetime study of excitations of the single giant system 
using the Born-Infeld action see [49]. The most important result from [49, 30] for us, is 
that all the deformations of the single threebrane giant graviton that we consider, are 
supersymmetric . 

The mixing of these operators with restricted Schur polynomials that have three 
columns (or more) is suppressed by a factor of order This factor arises from the 
normalization of the restricted Schur polynomials: the three column restricted Schur 
polynomials (with one short column) have a two point function which is smaller than 
the two column restricted Schur polynomials by a factor of order ;^[30]. Thus, at 
large N we can focus on the two column restricted Schur polynomials, which is a huge 
simplification. The analog of the statement that for operators with a dimension of 
0(1), different trace structures do not mix is: at large N restricted Schur polynomials 
Xr,(t,s) with R a Young diagram with n columns, each of which has length of 0{N), 
do not mix with operators XR',{r',s') that have n' ^ n columns. The fact that the two 
column restricted Schur polynomials are a decoupled sector at large N is to be expected. 
Indeed, at large N these operators correspond to a well defined stable semi-classical 
object in spacetime (the two giant system). We expect that n column restricted Schur 
polynomials are also a decoupled sector at large N for the same reason. 
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4. Simple Projectors 



When all Young diagrams labeling the restricted Schur polynomial have at most 2 
columns, the projector PR-^(^r,s) simplifies dramatically. The goal of this section is to 
explain the simplification and exploit it to efficiently build the relevant projectors. 

The projector PR-^(r,s) projects from representation R of Sn+m to representation 
(r, s) of Sn X Sm- One issue, which comphcates things considerably, is that repre- 
sentation (r, s) can be subduced more than once when irreducible representation R is 
decomposed into irreducible representations of the Sn x Sm subgroup. Giving a general 
rule to specify precisely how this multiplicity is resolved is nontrivial. For the operators 
we consider in this article, this problem does not arise. As soon as a system of three or 
more gaint gravitons are considered, it will be necessary to deal with this issue. For- 
tunately a well defined approach to resolving these multiplicities has been outlined in 
[10]. Basically, [10] considers elements in the group algebra CSn+m which are invariant 
under conjugation by CSn x CS„i- The Cartan subalgebra of these elements are the 
natural generalization of the Jucys-Murphy elements which define a Cartan subalge- 
bra for S'„[50]. The multiplicities will be labelled by the eigenvalues of this Cartan 
subalgebra [10]. It would be very interesting to work out the details of this proposal in 
the context of multigiant systems. 

Given R we imagine removing some boxes; after removing these boxes one is left 
with r. The removed boxes are assembled to produce s. If we specify both R and the 
boxes that are to be removed to obtain r we obtain every representation exactly once. 
For example, consider all possible operators that can be constructed from the following 
R by removing the three boxes shown 













* 




* 


* 





For the boxes that are removed, we must respect edges that are joined, which means 
the two boxes removed from the short column must remain stacked on top of each 
other. Thus, there are two possible irreducible representations s that can be produced, 
implying we can subduce two possible irreducible representations of SgxSs, by removing 
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these three boxes 



Thus, in total we'd get the following 5*9 x 5*3 irreducible representations subduced from 
the R given above (sum over all possible ways to remove boxes to get this result) 



The dimensions of these representations are 42, 48, 96, 27, 54 and 8 respectively. This 
sums to give 275 which is indeed the dimension of R. 

To construct the actual projector, we need only build an operator which will as- 
semble the removed boxes in the correct way to produce s. We will give an example 
of how to construct this operator; the general case should be clear. Lets start with the 
representation R shown, removing the boxes indicated 



Our projector will act in the subspace spanned by the two sets of states 



|i) = 



iT 



|2) 



We are using a Young- Yamonouchi basis. Thus, each state given above could be any one 
of = 3 states, corresponding to the number of ways to complete the labels. When 

acting in the subspace, the operator which organizes the boxes into representation s is 



ml 



J2 Xs{(^)TR{a) 



aeSr, 
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All that remains is to supply a formula for the action of Tji{a) (for a = 1, (12)) when 
acting on the subspace spanned by |1) and |2). The action of I^Ricr) on any Young- 
Yamonouchi state is well known. For the states above 

r^((i2))|i) = -i|i) + ^|2) 

r«((12))|2) = i|2) + ^|l) 

so that 

((12)) = -i|l)(l| + ^|2)(1| + i|2)(2| + ^|1)(2| 

and 

r«(i) = |i)(i| + |2)(2|. 
5. The Radial Direction 

In this section we describe a limit in which the dilatation operator simplifies signif- 
icantly. There are two columns in the Young diagrams labeling the restricted Schur 
polynomials. When the first column contains 0{y/N) boxes more than the second, 
the dilatation operator simplifies to a lattice realization of the second derivative. The 
Young diagram label itself defines the lattice. 



5.1 Three Impurities 

The three impurity operators arc built using many Zs and three Vs. To specify these 
operators, wc need to give the three Young diagrams labeling the restricted Schur poly- 
nomial. The second Young diagram, r, (which specifics an irreducible representation of 
Sn) is specified by stating the number of rows with two boxes (= 6o) and the number 
of rows with a single box (= bi). The third Young diagram label, s, (which specifies an 
irreducible representation of S3) and the first Young diagram label, R, (which specifies 
an irreducible representation of Sn+3) can now be built from r by specifying which 
boxes in R are to be removed to obtain r and how these boxes are to be organized into 
an S3 irreducible representation. There are 6 possibilities 



XAibo, h) = xvn {Z, Y) xb(&o, h) = xm (Z, Y) 
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Xcibo, bi) = xm (Z, y) Xnibo, 61) = Xffl_ (Z, Y) 



P 

XE{bo, bi) = xm {Z, Y) xf(&o, &i) = Xmm {Z, Y) 



The corresponding normalized operators are denoted using the capital letter O. In view 
of the discussion of section 3, we know that 60 is 0{N) and 61 ranges from or 1 to 
0{N). The action of the dilatation operator is given in Appendix A. 

The 7^-charge of an operator in the field theory maps into the angular momentum 
of the dual string theory state. Thanks to the Myers effect [51] the angular momentum 
of the string theory state determines its size. Identifying the two columns of the Young 
diagrams with the two threebranes, the number of boxes in each column determines the 
angular momentum and hence the size of each threebrane. In the limit that N — bo — 
0(N), bo — 0{N) and bi — 0{-\fN) we have non- maximal giants which are separated 
by a distance of 0{1) in string units. In this limit, we expect the dynamics to simplify. 
Indeed, the system should be described by two D2> brane giant gravitons with open 
strings stretching between then. The action of the dilatation operator becomes 



DOA(boM)^9YM {N-bo)xO[- 



DOeibo, bi) = --qIm {N - bo) [Ob(6o + 1, 61 - 2) - 2 Ob(6o, &i) + Osibo - 1, 61 + 2)] 
2\/2 

+^9ym (N - bo) [Oc(bo + 1, 61 - 2) - 20c{bo + 1, 61) + Oc(&o - 1, &i + 2)] 



2\/2 

DOcibo, 61) = -^9^m{N - bo) [Osibo + 1, 61 - 2) - 20b(&o, &i) + Ob(6o - 1, &i + 2)] 
—gYMiN - bo) [Ocibo + 1, 61 - 2) - 20c{bo, 61) + Oc(6o -l,bi + 2)] 



DOoibo, bi) = --gl-MiN - bo) [Ooibo + 1, 61 - 2) - 20b(&o, &i) + Od(6o -l,bi + 2)] 
2V2 2 



+^9YMiN - bo) [OE{bo + 1, 61 - 2) - 20^(60, 61) + 0^(60 -l,bi + 2)] 
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DOE{bo, h) = -^gl^uiN - bo) [Od(&o + 1, 6i - 2) - 20^,(60, &i) + Ooibo - 1, 61 + 2)] 
-IdYMiN - bo) pEibo + 1, 61 - 2) - 20e(6o, 61) + Oe(6o - 1, &i + 2)] 



i?OF(&o, &i) = Qym {N -bo)xO 



1 

A, 

These results have a natural interpretation. Notice that there arc four operators for 
which the 5"^ representation is the totally antisymmetric representation. We will see 
that there are also four operators for which the corresponding states remain super- 
symmetric; this agreement between the number of operators for which the Sm rep- 
resentation is the totally antisymmetric representation and the number of supersym- 
metric states, holds in general for the two giant system. Looking at the labels, it is 
natural to interpret 0^(6o,&i) as a state in which we deform only the larger three- 
brane. Recall from section 3, that deforming a single threebrane gives us a super- 
symmetric state so it seems natural for 0^(6o,&i) to remain supersymmetric. Sim- 
ilarly, Opibojbi) can be interpreted as a state in which we deform only the smaller 
threebrane and a similar comment can be made. The fact that the combinations 
Osibo^bi) + \^Oc{bo,bi) and Onibo^bi) + \^OE{bo,bi) are annihilated by D im- 
plies that there are another two supersymmetric ways to deform the pair of three- 
brancs. Finally, notice that if wc set 0^(60,^*1) — Oc{bo,bi)/y/2 = OB-c{bo,bi) and 
Ooibo, bi) - Ob(6o, bi)/V2 = Oo-Eipo, 61) we have 

DOB-c{bo, bi) = -2gl-M{N - bo) [OB-c{bo + 1, 61 - 2) - 20b_c(6o, 61) + Ob-c(&o -!,&! + 2)] 
DOD-E{bo: &i) = -'2gYM{N - bo) [OD-E{bo + 1, 61 - 2) - 20DMbo: &i) + OoMbo - 1, &i + 2)] 

The right hand side again is a discretization of the second derivative. It is the Young 
diagram itself that is defining the lattice. After recalling that the number of boxes in 
each column sets the angular momentum and hence the radius^ of the corresponding 
threebrane, its clear that the radius of the giant graviton together with local physics 
in this radial direction has emerged. 

5.2 Four Impurities 

For the case of four impurities there are nine possible operators that we can define 

XAibo, 61) = xm {Z, Y) XBibo, b,) = xm (Z, Y) 



''The giant graviton threebrane wraps an S"^ of a given radius. It is the radius of this S'' that we 
call the "radius of the threebrane" . 
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Xc{bo,bi)=xm iZ,y) XD{bo,h) = xm iz,Y) 



XE{hoM) = X 



{Z,Y) XF{bo,h) = xrn {Z,Y) 



Xcibo, 61) = xm_ {Z, Y) xnibo, 61) = Xm_ {Z, Y) 



Xi{bo,bi) = x^m iZ,Y) 



Again, the corresponding normalized operators are denoted using the capital letter O. 
The action of the dilatation operator is given in Appendix A. 

In the limit that N — bo = 0{N), bo = 0{N) and 61 = 0{\/N) the dynamics again 
simphfies. The action of the dilatation operator becomes 



DOA{bo,b^)^{N-bo)glMxO{- 



bi 



DOsiboM)^ 



—gYM (N - bo) [OB{bo + 1, 61 - 2) - 2 Ob(6o, &i) + Ob(6o - 1, &i + 2)] 
+-^9ym {N - bo) [Oc{bo + 1, 61 - 2) - 20c{bo + 1, 61) + Oc(6o - 1, 61 + 2)] 



DOcibo, 61) = -Y9YMiN - bo) [Osibo + 1, 61 - 2) - 20b(&o, &i) + Ob(6o - 1, &i + 2)] 

-IgYMiN - bo) [Oc{bo + 1, 61 - 2) - 20c(6o, 61) + Oc(6o - 1, 61 + 2)] 



DOD{bo,b,) 



-2g^^ {N - bo) [Ooibo + 1, 61 - 2) - 2 Ooibo, &i) + Oz5(&o - 1, &i + 2)] 
+-^^2 ^ (AT - bo) [OE{bo + 1, 61 - 2) - 20^(60 + 1, 61) + 0^(60 - 1, 61 + 2)] 
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DOE{bo, bi) = -2g'y^{N - bo) [Os(6o + 1, 61 - 2) - 20i^(6o, 61) + Oij(6o - 1, 61 + 2)] 

+ ^g^^{N - bo) [Ooibo + 1, 61 - 2) - 20^(60, 61) + 0^(60 - 1, 61 + 2)] 

+ '^9YMiN - bo) [Opibo + 1, 61 - 2) - 20^(60, 61) + 0^(60 - 1, 61 + 2)] 



^0^(60, 61) = -2glM (N - bo) [Op{bo + 1, 61 - 2) - 2 0^.(60, b,) + Op{bo - 1, 61 + 2)] 

+ '^9ym (N - bo) [OE{bo + 1, 61 - 2) - 20s(6o + 1, 61) + Os(&o -!,&! + 2)] 



DOcibo, bi) = -\glu (N - bo) [Og(6o + 1, 61 - 2) - 2 Og(6o, b,) + Og(6o - 1, &i + 2)] 

+^9ym (N - bo) [Onibo + 1, 61 - 2) - 20j,(6o + 1, &i) + Oi/(&o - 1, &i + 2)] 

i^O//(feo, &i) = -^^yM - ^0) [Onibo + 1, 61 - 2) - 2 Oj,(6o, 61) + Onibo - 1, 61 + 2)] 
/3 

+^91m (N - bo) [Ocibo + 1, 61 - 2) - 20g(&o + 1, bi) + Og(6o - 1, 61 + 2)] 

2 



DOj{bo,b,) = {N-bo)gi.MxO 



1 

A, 

We can again identify combinations of operators that are annihilated by D, that is, 
that are BPS. Apart from Oyi(&0!^i) and 0/(6o,&i) we have 0^(605 &i) + V^Oc{bo,bi), 
OoiboM) + V30s(6o,6i) + 720^(60,^1) and OaipoM) + y30H(&o,&i)- Notice that 
all of the BPS operators from this section and the last can be written as 

0-Bvs{R,r) = \fdsOR,{r,s){boM) 

s 

where ds is the dimension of the irreducible representation s of the symmetric group. 
Finally, notice that if we set -\/3(9s(6o, ^i) — Oc{bo, &i) = Os-cibo, bi), V20D{bo, &i) — 
0^(60,61) = OD-F{boM). OD{boM) - y/SOE(bo,bi) + V20F(bo,bi) = ODF-E{boM) 
and v^Og(&o, ^'i) - Onibo, &i) = Oo-Hibo, bi), we have 

DOB-c{bo, bi) = -2g^M{N - bo) [OB-c{bo + 1, 61 - 2) - 20b-c(6o, 61) + Ob-c(6o - 1, 61 + 2)] 

DOD-F{bo, hi) = -2g^M{N - bo) [Od_f(6o + 1, 61 - 2) - 2C»d_^(6o, 61) + OD-ribo -l,bi + 2)] 

DODF-E{bo, bi) = -'ig^MiN - bo) pDF-Eiho + 1, ^1 - 2) - 20DF-E{ho, bi) + ODF-E{ho -IM+ 2)] 

DOG-H{ho, bi) = -2g^j^{N - bo) [C»G-f/(^>o + 1, fei - 2) - 20G-H{bo, bi) + Oa-Hipo - 1, 61 + 2)] 

The right hand side is again a discretization of the second derivative. 
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6. Numerical Results 



In this section we describe the result of numerically diagonalizing the dilatation oper- 
ator. 

When setting up a numerical computation of the spectrum of the anomalous di- 
mension operator, we need to specify the maximum value for the difference between 
the number of boxes in the long column and the number of boxes in the short col- 
umn. Given this value, denoted amax; we are able to determine how many operators 
participate in our problem and wc arc able to describe the resulting spectrum rather 
explicitly. We will focus on the case of even Omax- In this case the difference between 
the number of boxes in the long column and the number of boxes in the short column 
is always an even number. 

6.1 Two Impurities 

For a given value of Omax there are 2 -|- 2ainax states in total. There are |amax + 1 zero 
eigenvalues (corresponding to supersymmetric states). The remaining eigenvalues are 

\ o 2 • • 1 o Omax ^ 

= SgYM^ « = 1, 2, • • • , + 1 . 

6.2 Three Impurities 

For a given value of Omax there are 1 -|- Somax states in total. There are 2ainax zero 
eigenvalues (corresponding to supersymmetric states). The remaining eigenvalues are 

\ o 2 • • 1 o Ctmax 

-^i = SfYM^ ^ = 1, 2, • • • , , 

each with a degeneracy of two and a single maximum eigenvalue A = 4aniax5'YM + 8fl'yM- 
This degeneracy almost certainly indicates a symmetry enhancement in the large N 
limit. 

6.3 Four Impurities 

For a given value of flmax there are 1 + |amax states in total. There — 1 zero 

eigenvalues. The eigenvalues are again evenly spaced with a level spacing of Sqym ^-^d 
they are again degenerate. The low lying eigenvalues 

\ o 2 • • 1 o Ctmax 

^'i = 8fyM« « = 1, 2, • • • , , 

have a degeneracy which alternates between 3 and 4. Thus, we find three eigenvalues 
A = 8gYM^ followed by four eigenvalues A = ISgyM: followed by three eigenvalues 
A = 24:gYM: followed by four eigenvalues A = ^"^gyM 
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If Omax is a multiple of 4, the larger eigenvalues are given by A = 4amax5'yM + 85'yM) 
A = ^a^^^g'^f^ + lQg'^f^ and 

Aj = 4:a^e^gYM + '^^9ym + '^^Wym « = 1, 2, • • • , ; 

all of these eigenvalues on the last line above arc non-degenerate. 

If Omax (which by assumption is even) is not a multiple of 4, the larger eigenvalues 
are given by A = 4amax5'yM + ^9ym with a degeneracy of 2 and 

Aj = 4amax5'Yjv^ + i^QYM + ^^Wym ^ = 1, 2, • • • , """^"^^ ; 

all of these eigenvalues on the last line above arc non-dcgcncratc. 

Once again, the degeneracies observed almost certainly indicate a symmetry en- 
hancement in the large N limit. 



7. Discussion 

In this article we have computed the one loop anomalous dimension of an operator 
built from 0{N) Zs and 3 or 4 y "impurities". What lessons can be learnt from these 
results, together with the results of [30], which dealt with the case of 2 impurities? 
Before we start the discussion, it is useful to recall the structure of the operators which 
participate in the case of two impurities 



XA{bo, bi) = xm {Z, Y) Xs(6o, bi) = xm 



iZ,Y) 



B 



XD{bo,bi) = xm 



[Z,Y) XE{bo,bi) = xmr^ iZ,Y) 



m BBB 



Firstly, the actual result for the one loop dilation operator (see Appendix A) is 
rather complicated. This is to be expected - it was obtained by summing a huge 
class of Feynman diagrams - much more than just the planar diagrams. On the other 
hand, the spectra of anomalous dimensions obtained are rather simple. To obtain a 
numerically tractable problem, we have been forced to keep the value of Omax finite 
(recall that this parameter measures the maximum value for the difference between the 
number of boxes in the long column and the number of boxes in the short column) . In 
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the large N limit^ Omax = go^s to infinity. Our discussion assumes we are working 
in this Omax oo hmit. 



We have restricted ourselves to a study of less than five impurities. This restric- 
tion is not imposed because our methods break down for five or more impurities, but 
only because the details of writing down the projection operators and evaluating the 
dilatation operator becomes increasingly complicated as the number of impurities is 
increased. In particular, there is no problem in principle with taking 0{N) impurities. 
The case of no impurities or one impurity are simple to handle analytically - all of these 
operators are annihilated by the one loop dilatation operator[30]. 

For the case of two impurities there are three times as many zero eigenvalue states 
as there positive eigenvalue states. There are positive eigenvalue states with a 
constant energy level spacing of Sqym- Thus, it is natural to associate an oscillator 
with a set of ~ states. With this assumption, the dilatation operator acting on the 
two impurity operators gives the spectrum of three harmonic oscillators with a level 
spacing of zero and a single harmonic oscillator with a level spacing of Sqym- Looking 
at the two impurity operators given above, we see that there are three operators with 
impurities in the antisymmetric representation (B) and one operator with the impurities 
in the symmetric representation (m). 

For the case of three impurities, the dilatation operator gives the spectrum of four 
harmonic oscillators with a level spacing of zero and two harmonic oscillators with a 
level spacing of ^gyM'^ ea.ch. oscillator again has ~ states. Looking at the three 
impurity operators given in section 5.1, we see that there are four operators (A,B,D 

and F) with impurities in the antisymmetric representation (§) and two operators (C 
and E) with the impurities in the W representation. 

For the case of four impurities, recall that we had an interesting degeneracy struc- 
ture - the degeneracy alternates between three degenerate states and four degenerate 
states. This is naturally explained as three oscillators with a level spacing of ^Oym 
a fourth with a level spacing of IQqym^ with this interpretation each oscillator again 
has ~ states. Thus, for the case of four impurities, the dilatation operator gives 
the spectrum of five harmonic oscillators with a level spacing of zero, three harmonic 
oscillator with a level spacing of ^Oym ^'^'^ with a spacing of IQqym- Looking 
at the four impurity operators given in section 5.2, we see that there are five opera- 
tors (A,B,D,G and I) with impurities in the antisymmetric representation (□), three 
operators (C,E and H) with the impurities in the §~' representation and one (F) with 



was defined in section 3. 
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impurities in the S representation^. 

It is rather easy to guess the result for a general number of impurities. If the number 
of impurities is even = 2n we expect to obtain a set of oscillators with frequency oui 
and degeneracy di given by 

uJi = 8igYM, di = 2{n-i) + l, i = 0,l,...,n. 

If the number of impurities is odd = 2n + 1 we expect to obtain a set of oscillators with 
frequency oui and degeneracy di given by 

uJi^SigYM, = 2(n - i + 1), i = 0, 

This conjecture passes a simple counting test: - di is equal to the number of restricted 
Schur polynomials that can be defined. Further, the degeneracies di match the number 
of each type of oscillator that can be defined: di is equal to the number of operators 
which have the impurities organized into a Young diagram with i boxes in the short 
column. 

A beautiful simple picture is emerging from the rather complicated formulas ob- 
tained for the dilatation operator: the dilatation operator is equivalent to a set of 
harmonic oscillators. For each type of operator there is a single oscillator and the 
frequency of the oscillator is determined by the representation which organizes the im- 
purities. Since a set of harmonic oscillators is an integrable system, this system we have 
studied here is an example of an integrable dilatation operator, obtained by summing 
planar and non-planar diagrams. 

J\f — A super Yang-Mills theory has an SU (4) 7?.-symmetry. A U (2) subgroup of the 
7?.-symmetry can be used to rotate the Ys and Zs into each other. Thus, although the 
operators wc have studied do not treat the Ys and Zs symmetrically, a more complete 
study working with the complete set of two column restricted Schur polynomials would 
yield a description in which Z and Y appear on an equal footing. If we are to recover 
this symmetric description its clear that the oscillators which emerge for each particular 
impurity configuration must themselves have frequencies which are multiples of Sqym 
- exactly as we have found. Exploring this further we should be able to show that 
the spectrum is filling out U (2) multiplets. By including more species of impurities we 
should also be able to see more of the expected global symmetry group. 

What is the AdS/CFT dual interpretation of our results? The operators we have 
considered are dual to giant gravitons. A connection between the geometry of giant 

^Note that by looking at the representation that organizes the impurities we have been able to read 
off the frequencies of the harmonic oscillators appearing. This is where it ends; in particular we are 
not claiming that Of are the operators corresponding to the frequency IQqym operators! Operators 
with a good scaling dimension are a complicated linear combination of the various possible Os. 



- 18 - 



gravitons and harmonic oscillators was already uncovered in [34, 35, 36]. This work 
quantizes the moduli space of Mikhailov's giant gravitons. Consequently one is cap- 
turing a huge space of states. It is this huge space of states that connects to harmonic 
oscillators. Our study has focused on a two giant system. Consequently, we know 
that the oscillators that we have captured are associated to this two giant system and 
excitations of it. We have thus arrived at a slightly more refined statement of how 
the harmonic oscillator enters. The critical reader might question whether our set of 
operators includes excitations corresponding (for example) to the two giant system plus 
a graviton. This would seem to be a small perturbation of the two giant system that is 
not an excitation of it - the graviton is an excitation of spacetime. We do not have such 
excitations among our states: these states correspond to operators with a small third 
column, which have decoupled at large N. In much the same way that by quantizing 
the possible excitation modes of a string one obtains a set of oscillators, its natural to 
think that our oscillators arise from the quantization of the possible excitation modes 
of a giant graviton. 

We have described a limit (the first column of the Young diagram contains 0{\/W) 
boxes more than the second) in which the dilatation operator simplifies to a lattice 
realization of the second derivative. It is the Young diagram itself that is defining 
the lattice. After recaUing that the number of boxes in each column sets the angular 
momentum and hence the radius of the corresponding threebrane, its clear that the 
radius of the giant graviton together with local physics in this radial direction has 
emerged. Notice that for BMN loops the number of lattice sites is 0{y/N); for the 
operators we have studied here the number of lattice sites is 0{N). 

One can contemplate a number of tests for our proposal. With a thorough un- 
derstanding of the vibrational modes of the giant graviton system, one could imagine 
realizing definite classical membrane geometries by considering coherent states of the 
oscillators. One might then compare the energy of these states with the volume of the 
membrane times the membrane tension. Although naively the field theory and gravity 
computations are never simultaneously valid, one might hope that for operators which 
are close to BPS, the perturbative result might agree with the strong coupling answer 
(see [39] for a relevant discussion). 

In addition to questions we pointed out above, there are a number of clear directions 
for further study. Given the simplicity of our results, it should be possible to construct 
an analytic solution. This is under investigation [52]. It would also be interesting to 
understand how our results are modified at higher loops. One could also consider the 
case of n > 2 column restricted Schur polynomials and more species of impurities. 
A much more general question would be to ask when (and how) simple systems are 
expected to emerge from multimatrix models. For a single matrix model, it is well 
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known that the planar hmit is nicely captured by the dynamics of N non-interacting 
non-relativistic fermions in an external potential. In this paper we have argued that the 
large N hmit of a class of operators dual to giant gravitons is captured by a collection 
of harmonic oscillators. Presumably every semiclassical object in spacetime (string, 
giant graviton, black hole, etc) is associated with the emergence of a simple system in 
the large limit of the corresponding class of operators in the field theory. Can we 
make this connection sharper and more useful? 
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Ives, Yusuke Kimura, Hai Lin, Dimitri Polyakov, Sanjaye Ramgoolam, Peter Roenne, 
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A. Dilatation Operator for Three or Four impurities 

The dilatation operator for the case of two impurities has been given in [30]. In what 
follows 
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B. A Discussion on Intertwiners 



When S"' acts on F®" n > 1 it furnishes a reducible representation. Imagine that this 
includes the irreducible representations R and T. Representing the action of cr as a 
matrix r{a), in a suitable basis we can write 
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If we restrict ourselves to an Sn-i subgroup of Sn, then in general, both R and S will 
subduce a number of representations. Assume for the sake of this discussion that R 
subduces R[ and R2 and that S subduces S[ and S2. Then, for a e Sn-i we have 
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Imagine that S[ — R[, that is, one of the irreducible representations subduced by R 
is also subduced by S. Then, a simple application of the fundamental orthogonality 
relation gives 
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where the form of the intertwiners has been spelled out. 
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